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We investigate theoretically the spin purity of single holes confined in vertically coupled 
GaAs/AlGaAs quantum dots (QDs) under longitudinal magnetic fields. A unique behavior is 
observed for triangular QDs, by which the spin is largely pure when the hole is in one of the dots, 
but it becomes strongly mixed when an electric field is used to drive it into molecular resonance. 
The spin admixture is due to the valence band spin-orbit interaction, which is greatly enhanced in 
Czh symmetry environments. The strong yet reversible electrical control of hole spin suggests that 
molecules with Gs-symmetry QDs, like those obtained with [111] growth, can outperform the usual 
G 2 -symmetry QDs obtained with [001] growth for the development of scalable qubit architectures. 


PACS numbers: 73.21.La,73.40.Gk,71.70.Ej 


Single spins confined in III-V semiconductor QDs are 
currently considered as potential qubits for solid-state 
quantum information processing, which combine fast 
optical and electrical manipulation with prospects of 
scalabilityji^— In the last years, heavy hole spin qubits 
have emerged as a robust and long-lived alternative to 
electron spins, as they can be less sensitive to dephasing 
from nuclear spins Significant advances have been re¬ 
ported on hole spin initialization, control and readout by 
means of optical excitations^TiiiiSr— , and different pro¬ 
posals for electrical control have been put forwardJ^^— 

Although most works so far have focused on QDs 
grown along [001], it has been noted that the C 2 point 
symmetry of such systems gives rise to a splitting of 
bright exciton states which limits the fidelity of optical 
hole spin preparation . No such splitting is however 
expected in [111] grown QDs owing to their higher (Ca) 
symmetryj22ii^ which hence become an alternative worth 
exploring. Early studies on single (In)GaAs/AlGaAs 
QDs grown along [111] have revealed that hole states 
have weak heavy hole-light hole (HH-LH) coupling due 
to the large aspect ratio, which is a prerequesite to ob¬ 
tain pure hole spins and minimize the impact of hyperfine 
interaction with the lattice nuclei;^ In turn, magneto¬ 
photoluminescence spectra have reported characteristic 
differences from [001] grown QDs which were ascribed to 
the influence of the C3 symmetry on the hole statesi^ 

In this paper, we move forward and study hole states 
confined in quantum dot molecules (QDM) formed by a 
pair of vertically stacked QDs grown along [111]. QDMs 
present several advantadges over single QDs for qubit 
development, including readout independency from ini¬ 
tialization and measurement protocols^, higher fidelity 
of spin preparatioiiiS and enhanced wavelength tun- 
ability with external electric fields, which greatly im¬ 
proves prospects of scalability/^ We consider [111] grown 
GaAs/AlGaAs QDMs with triangular shape, similar to 
those reported in Refs/ ^^i^^ , adding longitudinal magnetic 
and electric fields to control the Zeeman splitting and 
charge localization. Our calculations show that the HH 
spin purity is high when the hole is confined in individ¬ 


ual QDs, but severe spin admixture takes place when the 
electric field is used to form molecular orbitals. The spin 
admixture follows from the formation of orbitals with 
approximate point group symmetry, which enables 
otherwise forbidden spin-orbit interactions (SOI). The 
symmetry-induced SOI does not mix nearby Zeeman sub- 
levels, but it couples bonding and antibonding molecular 
states split by the tunneling energy. This is in sharp con¬ 
trast with usual [001] grown QDMs, with C 2 h symmetry, 
where tunneling is normally a spin-preserving process/^ 

Since the activation of SOI mechanisms is gener¬ 
ally associated with a descent of the system symmetry 
(e.g. system and bulk inversion asymmetry^^, QDM 
misalignment^), the enhancement of SOI for QDs 
is apparently counterintuitive. Yet, we observe strong 
spin admixture between hole states over 1 meV apart, 
which is 2.5 times greater than the largest spin-orbit an¬ 
ticrossing measured in [001] grown QDMs/^ We provide 
an explanation through group theory analysis of the the 
multi-band k-p Hamiltonian for holes, showing this is an 
exclusive property of C3 systems, and discuss the impli¬ 
cations of these findings for the development of hole spin 
qubit architectures. 

The Hamiltonian we use to describe hole states reads: 

H = Hsi? -|- Hb + M.strain + (Y(r) -|- 6 {(j)pz{v) — F z j) I. 

( 1 ) 

Here Mbf is the four-band Hurt-Foreman Hamiltonian^! 
for [111] grown zinc-blende crystals, which considers HH- 
LH subband coupling as in the Luttinger model^ but 
including position-dependent effective masses. Hs rep¬ 
resents the terms coming from a magnetic field applied 
along the growth (z) direction. Mstrain is the strain 
Hamiltonian, F(r) the band-offset potential, e the hole 
charge, (jjpz the piezoelectric potential, F an axial elec¬ 
tric field and I a rank-4 identity matrix. Further details 
on the Hamiltonian can be found in the Supplemental 
Material.— Hamiltonian CD is solved numerically after 
obtaining the strain tensors and piezoelectric fields us¬ 
ing the Comsol package. The eigenstates are Luttinger 
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spinors of the form: 

3/2 

N= E /I(r)|J = 3/2,J,), (2) 

Jz=-3/2 

where /" (r) is the envelope function associated to | J = 
3/2, J^), the periodic function with Bloch angular mo¬ 
mentum Jz. Jz = ±3/2 correspond to spin up and 
spin down HH components, while = ±1/2 corre¬ 
spond to LH components. The expectation value (Jz) = 

can be taken as a measure of the hole 
spin purity, with (Jz) ss ±3/2 indicating nearly pure HH 
spin up (fi) or spin down (JJ.) states. 

For our calculations we consider pyramidal 
GaAs/Alo.sGao.rAs QDMs with triangular confine¬ 
ment, similar to those obtained by metallorganic vapor 
depositionj^i^ The vertically stacked QDs are separated 
by a barrier of thickness d = 2 nm, although they are 
interconnected by a thin Alg.osGao.gsAs vertical quan¬ 
tum wire which enhances tunnel coupling -see structure 
in Fig. [TJa)-i^ A weak magnetic field oi B — 0.2 T is 
applied along the coupling direction. Figure [TJa) shows 
the energy of the first four hole states under the influence 
of a vertical electric field F. At F = 11 kV/cm we see 
two Zeeman-split doublets. The upper one (states |1) 
and 12)) corresponds to the main component of the hole 
spinor in the top QD, which is slightly bigger than the 
bottom QD. In turn, the lower doublet (states |3) and 
|4)) has the main component in the bottom dot. Since 
the increasing electric field favors the occupation of the 
bottom dot, a charge transfer anticrossing takes place 
at F = 14.7 kV/cm (resonant electric field, F^.), where 
bonding and antibonding molecular orbitals are formed. 
The behavior closely resembles that of [001] grown 
QDMs^^i^i^ except for one anomaly: the Zeeman 
splitting of both doublets is quenched near the resonant 
electric field, see Fig. (TJa) inset. 

The Zeeman splitting suppression can be seen as a van¬ 
ishing effective g-factor. Unlike in previous reports, how¬ 
ever, the origin of this effect cannot be ascribed to the 
different g-factor of the QD and barrier materials^ as 
the QDs and the vertical wire connecting them have sim¬ 
ilar composition. Because electrically tunable g-factors 
are of interest for spin manipulation^i^ we further inves¬ 
tigate into the origin of this phenomenon. Fig.jTJb) shows 
the hole Bloch angular momentum expectation value of 
the four states under consideration. As can be seen, away 
from the resonant field, (Jz) gradually approaches ±3/2, 
indicating that the hole states confined in individual QDs 
are nearly HH states with fairly pure spin. In the vicinity 
of Fr, however, (Jz) ~ 0, which means that for molecular 
states the spin becomes completely mixed. 

For comparison, since [111] grown QDs normally have 
either triangular or hexagonal shapej^i^li^ in Fig. [IJc- 
d) we study a QDM formed by hexagonal QDs. One can 
see that no spin mixing takes place near in this case. 
Actually, the behavior is now the same as observed in 
vertically aligned [001] grown QDMs with C 2 h symmetry 
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FIG. 1. (a) and (b): hole energy levels and Bloch angu¬ 

lar momentum expectation value of a triangular QDM grown 
along [111], as a function of a vertical electric field, (c) and 
(d): same but for a hexagonal QDM. Note the strong spin 
mixing for triangular QDMs near the resonant field Fr- The 
structures in (a) and (c) show the hole localization for each 
doublet. The inset in (a) shows the Zeeman splitting of the 
upper (solid line) and lower (dashed line) doublet. 


QDs, where tunneling is a spin preserving process;^ It 
follows that the strong spin mixing in Fig. [I](b) is not a 
consequence of the [III] crystal orientation, but rather 
of the triangular envelope confinement. In fact, we also 
observe it for triangles grown along [001], see Supple¬ 
mentary Material)^ It does not result from the presence 
or absence of strain either, as similar results are obtained 
using lattice mismatched materials such as InAs/GaAsi^ 
Likewise, it is not induced by the magnetic field, as Hs 
is a diagonal term which does not couple different spinor 
component!^ It does not take place in single QDs either. 
It is an exclusive property of QDMs with triangular con¬ 
finement. 

Further insight into the hole spin mixing mechanism 
is obtained in Fig. [51 which plots the weight of the four 
spinor components corresponding to each of the states 
jl) to |4) of the triangular QDM. Two conclusions can 
be extracted: (i) LHs play a minor role in all cases, the 
mixing is essentially between HH components with or¬ 
thogonal spin projections (Jz = ±3/2); (ii) states Jl) 
and j4) (panels (a) and (d)) seem to exhibit complemen¬ 
tary behavior, and so do 12) and 13) (panels (b) and (d)). 
This suggests that the spin mixing is due to independent 
interactions between jl) and j4) on the one hand, and j2) 
and j3) on the other. 

To understand the origin of the HH spin mix¬ 
ing we resort to a point group theory analysis. 
Eq. II]) Hamiltonian can be simplified as H « 
± ± (f^(r) ±eFz)) I, where we have disre¬ 

garded strain terms -which are weak for GaAs/AlAs 
heterostructures- and approximated Hsi? by the (con- 
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FIG. 2. Weight of the different Luttinger spinor components 
of the hole states in Fig.[TJa), as a function of the electric field. 
Panels (a) to (d) correspond to states |1) to |4), respectively. 
The states are almost exclusively HH (Jz = ±3/2). 


slant mass) Luttinger-Kohn Hamiltonian: 
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( 3 ) 


with 

P±Q = (71 ±73)(fc2 ±fc2) ;p 273)fc2 

R = + 273)^- + ^(72 - Ipk + k ^ ( 4 ) 

S = -^(72- 73 )fc+±(272 ± 73 )A:-fcz 


where 71 , 72,73 are the Luttinger parameters and k± = 
kj:±iky. One can then see that H has C 3 point symmetry 
set by the confining potential V{r). Near the resonant 
electric field, however, an additional approximate sym¬ 
metry must be considered. Even if the two QDs forming 
the QDM are not identical, the electric field restores an 
effective parity symmetry, the bonding and antibonding 
HH molecular orbitals forming even and odd functions 
with respect to a mirror plane in between the QDsj^ 
The corresponding point group is then Cah- Note that for 
to hold exact C^h symmetry, we need to impose the 
so-called axial approximation (72 = 73 ), which is actually 
valid for many HI-V materials, such as GaAs. Actually, 
we do not impose exact symmetry in the numerical cal¬ 
culations and nevertheless, the obtained results reveal, as 
expected, a high degree of symmetry. 

The anticrossing of Fig. [IJa) can be rationalized con¬ 
sidering the symmetry of the hole spinors in the double 
group C^h- Within this group |1) and |4) have E _^/2 
symmetry, while |2) and |3) have A 3/2 symmetry)^ The 
different symmetry of |1) and |2) (|3) and |4)) explains 
the lack of interaction within the Zeeman doublets, in¬ 
spite of the quasi-degeneracy. It also becomes clear why 
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FIG. 3. (a) Diagram of single band hole energy levels under 
a longitudinal magnetic field for a QDM with symmetry C 35 ,. 
The labels (A^, Jz) indicate the symmetry of the levels. X 
represents the rotational symmetry of the envelope function 
(A, E±) while p = 1 or n represents the even/odd parity and 
Jz indicates the non-zero component of the spinor. Thick 
arrows denote the symmetry allowed level couplings for the 
ground state, (b) Energy structure of anticrossing hole states 
before and after switching on off-diagonal terms in '. (c) 
Typical dissociation energy spectrum of holes in QDMs at 
B = 0 T, showing a bonding-antibonding ground state rever¬ 
sal at dc. (d) Calculated spin purity of states |1) and |2) for 
different interdot distances at the resonant electric field. Spin 
mixing is strongest around dc. 


|1) and |4) (|2) and |3)) interact separately, as observed 
in Fig. 

The above picture differs from the widely studied [001] 
grown QDMs with circular confinement, where the sym¬ 
metry of all four spinorial states involved in the molecular 
anticrossing is different^!, which results in the absence 
of interaction and hence spin-preserving tunneling. Siz¬ 
able hole spin mixing has been observed only in QDMs 
with significant misalignment's, because the symmetry is 
then completely reduced {Ci point group). Nevertheless, 
the largest spin anticrossing measured for such system is 
0.4 meV, corresponding to InAs QDMs with anomalously 
large lateral offsetjSS This is 2.5 times smaller than the 1 
meV gap we estimate in Fig. [1] and 6 times smaller than 
the 2.5 meV we predict for InAs/GaAs QDMsi^ 

To explain the unusual strength of the spin-orbit cou¬ 
pling in triangular QDMs, we can examine the envelope 
symmetry of (within axial approximation) and the 

ensuing eigenfuctions. In the C^h point group, the terms 
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of Eq. ([3|) form basis of the following irreducible repre¬ 
sentations: 


/A' E'!_ E'j^ Q \ 

E'l A' Q E'^ \ 

E'_ 0 A' E'l \ 

\0 E'_ E'l A' j 


( 5 ) 


A remarkable consequence is that for any hole state |n), 
the envelope functions of the spin up and down HH com¬ 
ponents, f- 3/2 ®’ must have the same 

symmetry except for the even/odd parity (e.g. A' and 
A"). This is a key factor in determining the strength 
of the spin admixture, as we show in the perturbative 
analysis below. 

Hamiltonian El can be split as a sum of diagonal and 
off-diagonal terms, H = Hq -fll', the latter being respon¬ 
sible for band coupling. If we disregard H', the levels 
anticrossing at E^ are those represented at the top of 
Fig. Ha), namely two Zeeman doublets formed by HHs 
with opposite spin (J^ = ±3/2) but the same rotational 
symmetry (A) and parity (' or "). Each doublet is split 
by a Zeeman term and separated from each other 
by an amount 2t, where t is the HH tunneling integral. 
Considering H' as a perturbative term, the mixing be¬ 
tween the spin up ground state \k^°'>) = (A',-1-3/2) and 
any spin down HH state is given by: 


i^k 





|*(°)) ( 6 ) 


where is the j-th intermediate state and Ej its 

corresponding energy. Notice that the strength of cou¬ 
pling is inversely proportional to AEhh = — E^, i.e. 

the energy difference between the spin up and down HH 
states. As explained with detail in the Supplemental 
Material,— the symmetry of H' operators, off-diagonal 
terms of Eq.([5|), translates into selection rules which 
make the numerator of Eq. (H vanish for all except the 
two paths plotted with thick arrows in Fig. Hs-)- Both 
paths involve excited LHs as intermediate states^, and 
the spin down HH is |j) = (A", —3/2), i.e. a state partic¬ 
ipating in the molecular anticrossing. This implies AEhh 
is small (few meV at most), and is in contrast with other 
point symmetries, where selection rules lead to coupling 
with higher excited HH states, so that AEhh is much 
larger. For example, if we consider QDMs with circu¬ 
lar QDs (point group Cook), the HH components coupled 
by H' no longer have the same rotational symmetry, but 
they differ by three quanta of azimuthal angular momen¬ 
tum Mj;— Consequently, there is no coupling between 
the Mz = 0 HHs forming the molecular anticrossing— 
Having a QDM structure is also essential, as then A! and 
A" are roughly split by the tunneling energy 2t, which 
can be made small enough for the SOI to be efficient. By 
contrast, in single QDs the strong vertical confinement 
would lead to several meV splitting. 

The suppression of the Zeeman splitting observed in 
Fig. Ha) can be also understood from the perturbative 


analysis. As indicated in Fig. Hb), the band coupling 
occurs between HH states belonging to different dou¬ 
blets. Because AEhh is smaller for the innermost states 
(2t — Ag) than for the outermost ones (2t -|- A^), the 
interaction is stronger, leading to an effectively reduced 
Zeeman splitting, Ag. 

It is clear from the discussions above that tunneling 
must be an important parameter to control the strength 
of the spin mixing. One might then expect that spin 
mixing is enhanced for long interdot distances d, when 
tunneling energy t is small. Fig. Hd) shows the spin 
purity of the ground state HH for QDMs with different 
d at resonant electric field. Interestingly, the maximum 
spin mixing is found at intermediate distances, d k. 2 nm. 
This follows from the characteristic, non-monotonous de¬ 
cay of hole tunneling in QDMs— As shown in the 
schematic of Fig. He), there is a critical distance dc where 
bonding and antibonding hole states are reversed. At 
this point, t has a relative minimum combined with large 
wave function delocalization, which enables the strong 
spin mixing. For d < d^ t increases rapidly, reducing 
the interaction. For d > dc, t eventually decreases but 
so does the wave function delocalization. As a result we 
gradually retrieve the single QD limit, were spin mixing 
is weak. 

Electrical control of hole spins in QDMs has been 
proposed as a key ingredient for scalable qubit 
architectures— So far, however, only [001] grown QDMs 
have been considered, where the main source of spin 
mixing was misaligment between the vertically stacked 
dotS)^ which is a difficult parameter to regulate exper¬ 
imentally. The Ca^i-symmetry-induced spin mixing de¬ 
scribed here arises as a more robust and manageable 
mechanism. It can also help increase the fidelity of spin 
control gate operations, as this requires the spin mixed 
states to (i) be energetically well resolved, and (ii) be 
able to form indirect excitons with large optical dipole 
strength— As for (i), we predict strong mixing between 
states 1 meV away from each other, larger than any 
previous measurement. As for (ii), unlike in misaligned 
QDMs, the spin mixing we describe is strong at the reso¬ 
nant electric field, where direct and indirect excitons have 
comparable optical strength. Another advantadge is the 
possibility to use weak magnetic fields, which limits the 
influence of the g-factor inhomogeneity of different QDs 
in the qubit scaling. 

In summary, we have shown that triangular QDs can 
be used to build QDMs with electrically controllable hole 
spin. The hole spin is well defined inside the individ¬ 
ual QDs, but the formation of delocalized molecular or¬ 
bitals with Csh symmetry enables SOI induced mixing 
with unprecedented strength. The reversible control, 
the strength of the interaction and the robust nature 
of the spin mixing mechanism imply that holes in tri¬ 
angular QDMs, like those obtained with [III] growth, 
form a promising system for quantum information pro¬ 
cessing with some advantadges as compared to circular 
[001] grown QDMs. 
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I. THEORETICAL MODEL AND PARAMETERS 


The Hamiltonian we use to describe hole states in QDs grown along the [001] direction 
is (atomic units): 


g[[CMi| ^ 5,(0011 Hgon ^ (^iooil(^) _ 5, 5 


JOOl] 


( 1 ) 


Here is the four-band Hurt-Foreman Hamiltonian^ for zinc-blende crystals, which 


considers HH-LH subband coupling and position-dependent effective masses: 

^ x.y.z ^ ^ ^ x.y.z ^ ^ 

u[001] _ ( 1 ^ ^ u u 

^BF ~ (9 9 kijlo 2 _^ 


+ ^{k:c{N - N')ky - ky{N - N')k,)I, 


1 


x,y,z 


Y,iki{N + N')kj + kj{N + N')ki){SiJj} 


i<j 


+-[{k+{N - N')k, -k,{N- A')fc+) J_ + {k,{N - N')k_ - k-{N - N')k,)S+] (2) 

with I a rank-4 identity matrix, {A, B} = \{AB + BA), X± = ± Xy and Jj the angular 

momentum i-component matrix: 
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By inserting 

eq. 

|3]in 

eq. 

0we end up 
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where, 


P' = + M)K + ky{L + M)ky + K{2M)K) + ^^{k^{N - N')ky - ky{N - N')k,) 

F” = lik^iL + 5M)k^ + ky{L + bM)ky + 2k^{2L + M)k,) + l{k^{N - N')ky - ky{N - N')k^) 

F'* = \{k^{L + bM)k^ + ky{L + 5M)ky + 2k,{2L + M)k,) - l{k^{N - N')ky - ky{N - N')k^) 

F* = \{k,{L + M)k, + ky{L + M)ky + k,{2M)k,) - i{k,{N - N')ky - ky{N - N')k,) 

R = ^[K{L - M)k, - ky{L - M)ky - t{k,{N + N')ky P ky{N + N') k,)] 

F_ = -^[k_Nk, + k,N'k_] 

5l = -j^[k,Nk+ + k+N'k,] 

S+ = -j=^[k+Nk, + k,N'k^] 

si = -j^[k,Nk_ + k.N'k,] 

C = -l{k,{N - N’)k_ - k_{N - N’)k,) 

Ct = -l{k+{N-N')k,-k,{N-N')k+) 

(5) 


with L, M, N, N' being the Stravinou-van Dalen mass parameters. 

By setting the parameters constant, Eq. (01) Hamiltonian turns into the Luttinger Kohn 


Hamiltonian: The Stravinou-van Dalen parameters are then related to the Luttinger param¬ 


eters 7i, 72, 73 as L — M = —872, 3F -|- M = —271 — 572, N — N' = 1 + ji — 272 — 873 and 


N + N' = - 873 . 

represents the terms coming from a magnetic held applied along the growth (z) 
direction, F:- 


Hs = - 



(x'^ + y^) + ^{xky 




72)1 72 J 




( 6 ) 


with I a rank-4 identity matrix, the angular momentum z-component matrix (J = 8/2), 
K = 4/8 the effective g-factor- and fiB the Bohr magneton. 

®^itraL is the strain Hamiltonian, formally identical to Eq. (0]) with the products 
kikj replaced by the strain tensor component Cjj. e and (ppz are the hole charge and 
strain-induced piezoelectric potential -which is also diagonal-, F an axial electric held and 
H(r) the band-ohset potential. For the triangular QDMs we use a structure like that 
obtained by metallorganic vapor deposition^. The structure is illustrated in Fig. [H using 
similar geometry and composition to the pyramidal QDMs of Refs.-*^. Thus, the QDs are 
made of GaAs, the barrier of Alo. 3 Gao, 7 As and the vertical wire connecting the dots of 
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FIG. 1: QDM formed by triangular QDs. Vertical electric and magnetic fields are applied. 

Alo.osGao.gsAs. We note that the vertical wire plays no critical role in the phenomena we 
describe. The robustness of the results are checked by carrying out additional calculations 
with triangular round edges and also breaking the symmetry with two round and a sharp 
edge (not shown). For comparison, we also study QDMs made of hexagonal QDs. The sides 
of the hexagons are taken to have the same dimensions as those in Fig. [H and a hexagonal 
wire is used which preserves the Gg symmetry. 


For [111] grown QDs, the hole Hamiltonian becomes: 

Hl‘“l = Hgp + + h 1‘“L + (l/(r) + e (4V‘l(r) - F 2 )) 1. (7) 

where we have obtained and from Eq. ([2]) by writing ki and J* appearing in 

the Hamiltonian as a function of the new coordinates k[ and J' according to. 


k = Mk' 


J = MJ' 


( 8 ) 


where M is the rotation matrix:- 


f 

V6 


1 

vT 
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_L J_ \ 

V2 Vs 
1 1 

y/2 Vs 

2 n 1 

3 0 V3 / 


(9) 


Once we reach the new matrices, the prime is removed from fcj for the sake of a better 
presentation. The Hurt-Foreman Hamiltonian in the [111] direction then reads: 
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where, 


P' = [Kill + 'l 3 )kx + fcy( 7 i + 'l 3 )ky + k^iji - 273)/c^] - i [kxi'yi - 272 - 873)- ^^(71 - 272 - 873) 

P” = [kxili - l 3 )kx + fcy(7i - l 3 )ky + fc;j(7i + 273)^2] - I [kxi'^1 - 272 - 873)fcy - /Cy(7i - 272 - 873)^3,] 

* 5 - = {(/ c _ 7 i / c 2 - k^jik-) + V 2 [kx{'y2 - l3)kx - ky{'^2 - J3)ky + i {kx{'^2 - l3)ky + ky{'-^2 - l3)kx)] 

-2 [2A;_(72 + 73)/c^ - k^'y^kJ\] 

S+ = {{k+jik^ - k;,jik+) + \pl [kx {'^2 - l 3 )kx - ky {'^2 - l 3 )ky - i {kx {'^2 - l 3 )ky + ^^(72 - l 3 )kx)] 

-2 [ 2 fc +(72 + 73)^2 - kz')^k+]] 

R = { k -{' l 2 + 273)^- - \/2 [^2(72 - 73)A:+ + /c+(72 - 73)^^]} 

C = -| [kz{'^i - 272 - 873)- fc_(7i - 272 - 873)^2] 

( 11 ) 

The Mb Hamiltonian rotation requires a special care. It was obtained in^ by initially 
disregarding the effect of the remote bands and later enclosing it, replacing the mass by 
the effective mass. Then, since we rotate the crystalline structure keeping the axes hxed 
and then the Bloch functions, the form of EIb does not change. However, the effect of the 
remote bands does change, so that one should replace 72 by 73 in the expression of the 
effectives masses. This modification must be introduced in to reach 

To calculate the strain eij in elastic constants of [001] grown het¬ 
erostructures, and rotate the axes. The resulting elastic constants are related 

to by: 

eg;;' = ^ (12) 

a,b,c^d 

Likewise, for the piezoelectric potential we rotate the axes and obtain: 

= (18) 

k 

with 

Hamiltonians m and ([7]) are solved numerically after obtaining the strain tensors and 
piezoelectric helds using the Comsol 4.2 package. The material parameters of GaAs, AlAs 
and InAs are taken from Ref.-, except for the crystal density, dielectric constant and piezo¬ 
electric coefficient, which are obtained from Ref.-. Linear interpolations are used for all 
alloys parameters. Luttinger parameters are inferred from the linearly interpolated masses. 
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II. CHARACTER TABLE AND PRODUCT TABLE FOR THE DOUBLE GROUP 


Cs/i 

The character table we use is: 
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J. 

K 

1 

u 

cn* 

1 

u 

cn* 

1 

u 

cn* 

1 

u 

cn* 

X + iy 

E'_ 

1 

u* 

u 

1 

u* 

OJ 

1 

u* 

u 

1 

u* 

u 

X — iy 

A” 

1 

1 

1 

-1 

-1 

-1 

1 

1 

1 

-1 

-1 

-1 

z 

E'i 

1 

u 

cn* 

-1 

—u 

-cn* 

1 

u 

w* 

-1 

—u 

-cn* 

Jx T iJy 

E'L 

1 

u* 

u 

-1 

-u* 

—OJ 

1 

u* 

u 

-1 

-u* 

—u 

Jx 

E-xji 

1 

—u 

u* 

i 

—ioj 

ioj* 

-1 

u 

-u* 

—i 

iu 

—iu* 

J-112 

Ei/2 

1 

-u* 

u 

—i 

iu* 

—iu 

-1 

u* 

—u 

i 

—iu* 

iu 

Jl/2 

E5/2 

1 

—u 

u* 

—i 

ioj 

—iu* 

-1 

u 

-u* 

i 

—iu 

iu* 


-E- 5/2 

1 

-u* 

u 

i 

—iu* 

iu 

-1 

u* 

—u 

—i 

iu* 

—iu 


E -312 

1 

-1 

1 

i 

—i 

i 

-1 

1 

-1 

—i 

i 

—i 


E3I2 

1 

-1 

1 

—i 

i 

—i 

-1 

1 

-1 

i 

—i 

i 



( 14 ) 


where ca = e* 3" and Ji is the i-th component of the angular momentum. 


The table of products of C^h irreducible representations is: 
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( 15 ) 


III. SYMMETRY OF THE HAMILTONIAN AND THE WAVE FUNCTIONS 


The employed Hamiltonian is mass-position-dependent. Actually, in our system 
the mass parameters are constant within every domain, and have a sudden jump at the 
edge between neighboring domains. For the sake of easiness, to discuss on symmetry, we 
consider the Luttinger-Kohn constant-mass parameters limit of EIff- 


The symmetry of the Hamiltonian (eqs. IU HU] employing constant mass parameters) 
including a triangular confining potential and an axial magnetic field is C^. However, within 
the axial approximatioiil^ (72 = 73) it reaches C^h- The symmetries of their matrix element 
operators can then be calculated from the above character table and the expressions on eqs. 
m [m (assuming constant mass parameters): 


^Hlk - 
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E'L 
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A' 
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E'_ 
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A' 

E'L 
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K 

A' 

/ 


( 16 ) 


Accordingly, the symmetry of the envelope bonding/anti-bonding ground state functions 
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must be: 


( (&) ^ 


( A" (a) \ 

EUa) 

and 

K(b) 

E'_ {b) 


E'L (a) 

V (a) ) 


\A'{b) ; 


( 17 ) 


where the labels (a), (6) indicates the bonding/anti-bonding character of the components. 
These envelope components, combine with the Bloch functions yielding the wave function. 


The Bloch functions are built as the symmetry-adapted product of J = 1 angular momentum 
functions and the J = 1/2 spin functions. However, the presence of the mirror symmetry 
ah allows to employ bonding ( xi^h) = 1) and anti-bonding ( xi^h) = —1) J = 1 angular 
momentum functions. For example, employing anti-bonding angular momentum functions 
we have: 


13/2,3/2) = -j.\(X + ;Y) t) |3/2,-3/2) = ^\{X --lY) 4 .) 

(18) 

|3/2,1/2) = t> -^\(X + tY) 4) |3/2, -1/2) = y^lZ 4) + ^|(X - i F) f) 

The bonding J = 1 angular momentum functions are like the antibonding where X, Y, and 
Z are replaced by Jx, Jy, and J^. As a result, the table of products allows us to determine 
the symmetries of the Bloch functions: 


(19) 


|3/2, 3/2“) —)• E^/2 |3/2,3/2^) —)• E_^j2 

|3/2, l/2“) —)■ E_^/2 |3/2,1/2'’) -)■ Eii2 

|3/2, -l/2“) ^ E^,2 |3/2, -1/2'’) ^ E_y2 

|3/2, -3/2“) ^ E_ 3/2 |3/2, 3/2'’) ^ E 3/2 

The Bloch functions symmetries required to combine with the envelope components eq. [T 7 ] 
(left) yielding i?3/2 and i?-3/2, and those required to combine with eq. [T 7 ] (right) yielding 
^3/2 and E_^i2, are: 


^ E 3/2 ^ 


^ E-312 ^ 

E-512 

and 

El 12 

E 5/2 


E-I/2 

V ^- 3/2 / 


V ^ 3/2 / 


^ E-3j2 ^ 


^ T^3/2 ^ 
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-E-5/2 

E-112 


E^I 2 

\ ^3/2 ) 


\ ^-3/2 / 


( 20 ) 
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IV. SPIN MIXING IN QDMS GROWN ALONG [001] 


In the paper we have considered triangular QDs grown along the [ 111 ] direction because 
they are formed naturally in that direction. The physics leading to hole spin mixing is 
however connected with the envelope symmetry, and does not depend on the crystal orien¬ 
tation. To illustrate this point, in Fig. [ 2 ] we plot the expectation value of the Bloch angular 
momentum (J^) as a function of the interdot distance d for the upper Zeeman doublet of 
GaAs/Alo.aGao.yAs QDMs (states | 1 ) and |2)). All the parameters are taken as in Fig. 3 (d) 
of the paper, except that now the QDM is grown along [ 001 ] -i.e. Hamiltonian instead 
of As can be seen, the picture is qualitatively the same as that obtained for [ 111 ] 



d (nm) 

FIG. 2: Calculated spin purity of states |1) and |2) for different interdot distances at the resonant 
electric field for a triangular QDM grown along [001]. 


grown QDMs (Fig. 3 (d) in the paper). The main difference is that the critical distance dc, 
where the bonding-antibonding reversal takes place and hole spin mixing is maximized, is 
now shifted towards longer interdot distances. This is because the effective masses of HH 
along [ 001 ] are lighter than those along [ 111 ]. Therefore, the tunneling is stronger. 

Indeed, if we take constant mass parameters in eq. lU we obtain the corresponding 

Luttinger-Kohn Hamiltonian: 
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0 ^ 
R 


( 21 ) 


= - 


f P + Q -S R 
-S^ P-Q 0 
R^ 0 P-Q S 
\0 R^ P + Q J 


with 


P + Q = [(71 ± l2){kl + kl) + (71 T 2-i2)kl)]/2 
R = [-\/372(A:^ -kl)+i 2 ^/ 3 'y 3 kxky ]/2 

S = 'y3'/^{k^-iky)k^ 

Also, if we take constant mass parameters in eq. [101 obtain: 


( 22 ) 


f P + Q -S 


LK 2 


R 


-5t 
R^ 

VO 


P-Q 0 

0 P- 

R^ 


0 \ 
R 
Q S 

P + QJ 


( 23 ) 


with 


F ± Q = (71 ± iP{kl + kl) + (71 T 2^,)kl 


R 

S 


= -^{72 + 2-f3)kl + ^(72 - -f3)k+k, 
- _72 


( 24 ) 


-^(72 - 73)^^ + ^ (272 + l 3 ) k - k . 

By comparing F + Q in Eq. ( 12 ^ with Eq. ( 12411 . one can note the different HH effective 
masses in the z direction: 1/(71 — 272) vs l/(7i — 273). 

It is worth noting that the R operator in does not have C3 rotational symmetry 

(the C3 character table can be obtained from that of C3/1 by just considering rotations and 
keeping rows 1 - 3 , 7-8 and 12). However, for many III-V materials including GaAs, 72 ~ 73 
and one can approximate them both by 7 = (72 + 73 )/2 in the R and R^ matrix elements, 
thus yielding with axial symmetry,— that is reduced to G3 (or C3h) symmetry by the 

triangular conhning potential and the magnetic field. On the other hand, Hamiltonian 
does have C3 symmetry, but -as discussed in the paper- the axial approximation in both R 
and S matrix elements is needed to display exact Csh symmetry. 

In short, in both crystallographic directions the Hamiltonian has approximate C3h sym¬ 
metry, which becomes exact if the axial approximation is assumed. Similar considerations 
on the axial approximation hold for the strain terms. 


10 






V. SPIN MIXING IN INAS/GAAS QDMS 


Next, we consider InAs/GaAs QDMs grown along [ 001 ], similar to those obtained by self- 
assembled growtbli but with triangular (pyramidal) QD shape. Unlike for GaAs/AlGaAs, 
no interdot wire is present in this case. On the other hand, strain and piezoelectricity now 
play a signihcant role. 

Figure [ 3 ] shows the spin purity of the four hst hole states as a function of an external 
electric held. One can see that also in this case there is a strong spin mixing (|(J^)| <C 
3 / 2 ). An inspection of the individual spinor components (not shown) reveals that also for 
InAs/GaAs LH components play a minor role. Most of the mixing in Fig. [ 3 ] comes from 
admixture between HH = - 1 - 3/2 and = — 3/2 components, following the symmetry- 
induced mechanism described in our paper. 



FIG. 3: (a) and (b): hole energy levels and Bloch angular momentum expectation value of a 
triangular InAs/GaAs QDM grown along [001], as a function of a vertical electric field. The insets 
in (a) show the hole localization for each doublet. 
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It is remarkable that the spin mixing here takes place for bonding and antibonding states 
split by more than 2 meV. Besides, the spin mixing takes place over a wide window of electric 
helds (wider than for GaAs). This is because of the stronger SOI of InAs as compared to 
GaAs.— These results suggest that the eventual design of triangular InAs/GaAs QDMs 
would also form a promising system for hole spin manipulation. 


VI. PERTURBATIONAL ESTIMATE OF THE SPIN MIXING STRENGTH 


In this section we expand the discussion of Fig. 3 based on perturbation theory. The goal 
is to show that the spin in triangular QDMs is much stronger than in circular QDMs. 

The Hamiltonian describing the hole states in a QDM, be it or can be split 

as: 

H = Ho + H'. (25) 


Here Hq are the diagonal terms, whose eigenfunctions are single-band HH or LH states: 




mH^ixn) = f-i/2{xn 


/i\ 

0 

0 

voy 

/o\ 

0 

1 

\V 


, = fi/2ixn 


, = f-3/2{xn 


/o\ 

1 
0 
voy 

/o\ 

0 


0 

VV 


( 26 ) 


where /j^(X^) is the envelope function of X rotational symmetry and p parity. In turn, 
H' represents the off-diagonal terms of H, coming from H^i? and H^jram- This term is 
responsible for the band coupling. Without loss of generality, because the symmetry of H^^;' 
and M.s-train is the same, in what follows we consider GaAs/AlGaAs QDMs, where Mstrain is 
negligible. The analysis is further simplihed replacing H^f by its constant mass analogue, 
Mlk, Eqs. or fl 23 |) within axial approximation. This leads to the following expression 
for H': 
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(27) 


H' 


^ 0 

-A 

R 

0^ 


0 

0 

R 


0 

0 

S 

1 0 



0/ 


Within this approximation the matrix element operators R and S are just proportional to 
and k_kz, respectively. 

Considering H' as a perturbative term, the mixing of states up to second order is given by: 


ihb = E (e 

iz^k \jz^k 




(0)\ /,Tf(0)| 


,(0)\ 


|4,f)) 


El - £)» 


Et - El 


(28) 


Note that the coupling between requires 

as intermediate state, yielding the contribution: 


(A'f‘|S|Jff){Xf|fl|X«) 

where AEhh is the energy difference between the HH states 
and AEih that between \'^HHii.i^k'°)) state 

Alternatively, this coupling can also be achieved with intermediate state, 

yielding the contribution: 


(A'f|fl|Af>(Af |S|A'f) 


(30) 


AEhhAEih 

The matrix elements in the numerator of Eqs. ([29]), fl30|) determine the selection rules in 
the band coupling process. E.g. in the C^h group, the matrix element operator S has E'!_ 
symmetry (i? has Ej,_), see Eq. (ITB]) . Then, a totally symmetric A' ground state 

must couple, via S', with a state of symmetry E" {A' ®E'!_ = A', otherwise the 

integral is zero). Next, with symmetry E" (the complex conjugate of that of 

must couple, via i?, with symmetry {E'^®E'^®A" = A', otherwise the integral 

is zero). Likewise, of A' symmetry (note that A' and A" are reals and then coincide 

with their complex conjugates) can couple via R with of E'_ symmetry. Then, 

of E'j^ symmetry will couple, via S, with of A" symmetry [E'j_ 0 E'!_ 0 A" = A'). 

The above reasonings lead us to define the C^h allowed couplings, represented by thick 
vertical lines on the left side of Fig. |H Blue and yellow arrows correspond to either con¬ 
tribution, eqs. fl29|) and fl30|) . In the hgure we have simplihed the notation of the states in 
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eq. (12^ to {XP,Jz), where X is the rotational symmetry of the envelope function {A, 
and in C 3 or in Cqo), p = / or // represents the even/odd parity, and Jz indicates the 
non-zero component of the four-fold spinor. 

One can now compare with the case of circular QDMs, where the group is C^oh and the 
envelope functions are labeled by and parity. In this group, the matrix element operator 
R (S) is even (odd) and has Mz = —2 (M^ = —1). Now the ground state of 

0' symmetry can couple via S (—1") with a of symmetry 1". Then, of 

symmetry —1" couple, via R (—2'), with of 3" symmetry (because —1 — 2 -|- 3 = 0 

and // $$ / 0 // = /). 

couple, via R (—2') with 2'. In turn, of —2' symmetry 

couples, via S (—!"); with of 3" symmetry. 

Taking into account Fig. [Hand Eqs. fl2^ and fl30|) we can see that AEhh in the denominator 
involved in C^h is quite smaller than that involved in Cook and, therefore, the interaction 
should be much stronger. 

As a matter of fact, the strong spin mixing at the resonant electric helds is a singular 
behavior of triangular QDMs. In Cnh symmetries with n > 3 the spinor fourth component 
{Jz = —3/2) has different rotational symmetry symmetry than the hrst one {Jz = 3/2), 
as in the above discussed case of Cooh- Then, for similar reasons, the coupling between 
the states belonging to the hrst two doublets are also forbidden. In the case of C 2/1 QDM 
the symmetries of S and R are Bg and Ag respectively. Bonding (anti-bonding) molecular 
orbitals are of Ag {B^) symmetry. Then, the doublet antibonding ground state (bonding 
hrst excited state) are B^ I/, Bu Ij- {Ag fl, Ag J/). Therefore, any coupling among these four 
states is forbidden by symmetry, as can be easily checked with the help of Eqs. fl2^ and 
(EOD. 
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with symmetry (left column) and Caoh (right column). The labels {XP,Jz) indicate the sym¬ 
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ground state, as inferred from Eqs. ([29]), (fSO)) . 

^ J. van Bree, A. Yu. Silov, P. M. Koenraad, M. E. Flatte, and C. E. Pryor, g factors and 
diamagnetic coefficients of electrons, holes, and excitons in InAs/InP quantum dots, Phys. Rev. 
B 85, 165323 (2012). 

^ Q. Zhu, K. F. Karlsson, M. Byszewski, A. Rudra, E. Pelucchi, Z. He, and E. Kapon, Hy¬ 
bridization of Electron and Hole States in Semiconductor Quantum-Dot Molecules, Small 5, 
329 (2009). 

® F. Michelini, M. A. Dupertuis, and E. Kapon, Novel artificial molecules: Optoelectronic prop- 


15 

































































erties of two quantum dots coupled by a quantum wire, in Proceedings of the Ifth International 
Workshop on Computational Electronics, Pisa, 2010, edited by G. Basso and M. Macucci, p. 
255. 

® Q. Zhu, J. D. Ganiere, Z. B. He, K. F. Karlsson, M. Byszewski, E. Pelucchi, A. Rudra, and 
E. Kapon, Pyramidal GaAs/AlzGal-zAs quantum wire/dot systems with controlled heterostruc¬ 
ture potential, Phys. Rev. B 82, 165315 (2010). 

^ J.B. Xia, Effective-mass theory for superlattices grown on (llN)-oriented substrates, Phys. Rev. 
B 43, 9856 (1991). 

® I Vurgaftman, JR Meyer, and LR Ram-Mohan, Band parameters for III-V compound semicon¬ 
ductors and their alloys, J. Appl. Phys. 89, 5815 (2001). 

® O. Madelung, Semiconductors: Data Handbook, (Springer-Verlag, Berlin, 2004). 

L.G Lew Yan Voon and M. Willatzen, The kp method: electronic properties of semiconductors, 
(Springer Science & Business Media, 2009). 

E.A Stinaff, M Scheibner, A.S Bracker, I.V Ponomarev, V.L Korenev, M.E Ware, M.E Doty, 
T.L Reinecke, and D Gammon, Optical signatures of coupled quantum dots. Science 311, 636 
(2006). 

The spin admixture we report is proportional to the extradiagonal elements effective mass, 72 
and 73 . The value of these parameters is in turn modulated by the strength of the spin-orbit 
interaction. Then, the extent of the spin admixture is ultimately due to spin-orbit interaction. 
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